Magnetoresistance and electronic structure of asymmetric GaAs/Alo.sGao.rAs double 
quantum wells in the in-plane/tilted magnetic field 
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Bilayer two-dimensional electron systems formed by a thin barrier in the GaAs buffer of a standard 
heterostructure were investigated by magnetotransport measurements. In magnetic fields oriented 
parallel to the electron layers, the magnetoresistance exhibits an oscillation associated with the de- 
population of the higher occupied subband and the field-induced transition into a decoupled bilayer. 
Shubnikov-de Haas oscillations in slightly tilted magnetic fields allow to reconstruct the evolution 
of the electron concentration in the individual subbands as a function of the in-plane magnetic 
field. The characteristics of the system derived experimentally are in quantitative agreement with 
numerical self-consistent-field calculations of the electronic structure. 

PACS numbers: 73.20.Dx, 73.40.-c, 73.50.-h 
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I. INTRODUCTION 

In an idealized, infinitely narrow two-dimensional (2D) 
system the in-plane component of the magnetic field cou- 
ples only to the electronic spin degree of freedom. In 
real samples, however, the finite size of a 2D system in 
the growth (z) direction often leads to strong orbital ef- 
fects of thfi_in-plane field. For example, 2D-2D tunnel- 
ing studieslil of a system consisting of two nearby nar- 
row quantum wells, which is the simplest structure with 
non-trivial growth direction degree of freedom, have re- 
vealed dramatic effects of the in-plane field on the Fermi 
surface topology. In this weakly-coupled double quan- 
tum well sample the in-plane field, -By, in effect displaces 
the origin of the two layer Fermi circles by |e|i3||(i/?i and 
allows a sweep of one past the other at a critical field 
Be ~ 2hkp/\e\d, where d is the layer separation and kp 
is the individual layer Fermi wavevector. 

The crossing of Fermi surfaces is replaced with more 
complicated patterns in samples with strongly coupled 
quantum wells in which the tunneling rate dominates 
the electronic scattering, rate within a well. In a sim- 
ple tight-binding modelcl a partial energy gap, equal to 
the bonding-antibonding gap at zero field, opens at the 
wavevector corresponding to the crossing point of the dis- 
placed dispersion curves for the uncoupled layers. While 
the upper subband above the energy gap maintains a 
nearly parabolic shape, a saddle point develops in the 
lower energy subband. With increasing the bottom 
of the upper subband moves above the Fermi energy at 
a critical field i?c,i, resulting in a sudden drop in the 
density of states. At the second critical field, i?c,2, the 
saddle point of the lower subband approaches the Fermi 
energy and the density of states diverges. The mag^ 
netoresistajice oscillation observed on coupled doubleO'cl 
and tripletl quantum wells represents a striking manifes- 



tation of the two distinct van Hove singularities in the 
-dependent density of states. For these samples the 
tight-binding model has provided an accurate quantita- 
tive estimate for the critical in-plane fields. Bilayer syi 
tems realized in wide single quantum wells also displa; 
the magnetoresistance oscillation, however, the orbital ef- 
fects of the in-plane magnetic field are more complex here 
and cannot be captured by the simple tight-binding ap- 
proximation. Instead, the nunicrical self-consistcnt-field 
technique has proven successfulj for wide quantum wells. 

The importance of orbital effects of the in-plane mag- 
netic field has recently been emphasizeclQ also in the coru 
text of metal- insulator transition studies in Si-MQSFETEl 
and GaAs/Alo.sGao.rAs single heterojuncions.Bil2l The 
observed dramatic response to iiL-both metallic and 
insulating phases has been attributedQ'liSl to the sniopeipu- 
pling to in-plane field as well as to the distortioiiEJil3 of 
carrier Fermi surfaces. A detailed quantitative under- 
standing of these phenomena requires to incorporate the 
orbital effects of -B|| using non-perturbative approaches. 

In this paper we present a study of magnetotransport 
properties of GaAs/Alo.sGao.rAs heterojunctions with 
an additional thia Alo.3Gao.7As barrier introduced into 
the GaAs buffer. The structures were specially designed 
to combine properties of single-junction and bilayer sys- 
tems. A comparison between theory and experiment in 
our samples constitutes an excellent test for the reliabil- 
ity of the numerical self-consistent-field technique applied 
to systems with complex growth direction geometries. In 
Section II we present measured resistance oscillation with 
in-plane field and Shubnikov-de Haas (SdH) oscillations 
recorded at magnetic fields slightly tilted from the 2D 
layer plane. Theoretical calculations of the B|| -dependent 
electronic structure are used, in Section III, for quanti- 
tative interpretation of the experimental data. Our con- 
clusions are summarized in Section IV. 
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II. EXPERIMENTAL RESULTS 

We studied modulation-doped GaAs/Alo.sGao.yAs 
heterojunctions grown by molecular-beam epitaxy on 
(100) oriented semi- insulating GaAs substrates and pat- 
terned in the standard Hall-bar samples. A thin barrier 
formed by eight monolayers of Alo.3Gao.7As was grown 
inside the GaAs buffer producing two coupled quantum 
wells, one of a distorted rectangular shape and the other 
one with a nearly triangular geometry, as shown in in- 
sets of Fig. |l|. We chose the thickness and position of the 
barrier to achieve high population of the second (anti- 
bonding) subband. The particular structure design was 
based on zero-field numerical self-consistent simulations 
that predict minimum bonding-antibonding gap for bar- 
rier positioned near the node of the first excited state 
wavefunction (peak of the lowest energy state wavefunc- 
tion) of the original single-junction. Note that unlike thCN* 
conventional bilayers realized in double quantum wellsO 
or wide single wells our system is intrinsicly, strongly 
anisotropic. At zero magnetic field the lower (bonding) 
subband wavefunction has dominant weight in the rect- 
angular well while electrons from the higher (antibond- 
ing) subband are more likely to occupy the triangular 
quantum well. Parameters of the two studied samples 
are summarized in Table ^ The partial occupations of 
the bonding and antibonding subbands at = were 
obtained by Fourier analysis of the low-field SdH oscilla- 
tions. The Hall measurement at low perpendicular mag- 
netic fields provided an independent check for the total 
carrier density and was used to determine sample mobil-vf 
ities. The insets in Fig. ^ show the confining potential^ 
profiles and wavefunctions of the occupied subbands for 
two samples. 



critical field i?c,i and the upper critical field Be, 2- The 
values of critical fields are i?c,i = 5.8 T, Bc.2 = 11.7 T for 
sample A and Bc,i = 5.9 T, Bc,2 = 10.7 T for sample B, 
respectively. The in-plane magnetic field dependence of 
subband occupations, Na{B\\) and iVf,(i3||), were obtained 
from SdH oscillation data measured at different field-tilt 
angles. Typical SdH traces are shown in Fig. ||. In gen- 
eral, the oscillations are not periodic in 1/B± {B± is 
the perpendicular component of the magnetic field) since 
both B± and vary when the field is swept at a fixed 
tilt angle. Therefore, the standard Fourier technique does 
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90 


2.34 


2.77 


0.76 


3.42 
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TABLE I. Sample parameters: w is the thickness of the 
rectangular well, A^;, and Na are the bonding and antibonding 
subband densities obtained from SdH measurement, Nuaii is 
the total 2D electron density derived from low-field Hall data, 
and A'^ is the average between the SdH and Hall values for the 
total density used in the numerical simulations. 

Magnetotransport data were collected at the temper- 
ature 0.45 K and in magnetic fields ranging from to 
23 T, using both dc and low-frequency (13 Hz) ac tech- 
niques. First, we measured the resistance R^x for a mag- 
netic field precisely parallel with the 2D layer plane. The 
recorded magnetoresitance, AR^x/Rxx — iRxx{B^^) — 
RxxiO))/RxxiO), is plotted in Fig. |]. Consistent with 
the bilayer nature of the studied samples we observe an 
oscillation on the magnetoresistance trace with the lower 
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FIG. 1. The magnetoresistance trace (solid line) recorded 
for samples A and B at in-plane magnetic fields. For both 
samples the magnetoresistance curve is multiplied by a con- 
stant to fit the calculated density of states (dashed line) at 
high By . The insets show the band-profile and the wavefunc- 
tions of the bonding and antibonding subbands at i? = 0. The 
wavefunctions are shifted in the vertical direction so that their 
asymptotic values match the corresponding quantum levels 
on the energy axis. In both samples, the levels lie below the 
Fermi energy, Ef, indicating two occupied subbands at zero 
field. 

not apply. Instead, we use the measured distance be- 
tween valleys surrpimding individual peaks to identify 
subband densities.tj For < i?c,i one type of oscilla- 
tions, corresponding to the lower density (antibonding) 
subband, can be detected at small tilt angles while two 
subbands are clearly visible at higher angles. With in- 
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creasing B\\ the distribution of electrons between the two 
subbands changes, as shown in Fig. ^. EventuaUy, the 
antibonding subband is depopulated at _B|| — Bd- At 
higher in-plane fields two regimes can be distinguished. 
For _Bc,i < < -Be, 2 the SdH oscillations are periodic S 
indicating a single occupied subband. Consistently, the 
corresponding density is equal to the total 2D density 
measured at zero in-plane magnetic field. The character 
of the oscillations changes abruptly at _B|| — Bc.2] the 
low in-plane field non-periodicity is recovered and the 
obtained carrier density is significantly reduced (see Fig. 

These results suggest that the Fermi sea splits into 
two disconnected parts of densities that vary with . 
Only one Fermi surface, however, can be detected from 
our SdH data. 

800 



600 - 



400 



200 



G 



Sample A 


B^,=5.8T 
















2.4r ^^^^^^^^^^ 









12 



16 



20 



24 



B„ (T) 



600 



400 - 



200 



Sample B 


B,=5.9T 




6.74° A A 






4.54° 






' 3.7f 



















12 



16 



20 



24 



B„ (T) 



FIG. 2. Typical magnetoresistance traces measured for 
small angles between the sample plane and the magnetic field 
direction. The curves for higher angles are shifted upward by 

75 n. 



III. THEORETICAL ANALYSIS 

The single-particle Hamiltonian for an electron con- 
fined in a GaAs/Alo.sGao.yAs heterostructure and sub- 
jected to in-plane magnetic field can be written as 



H 



1 



2m 



-{p+\e\Ay + Vconfiz), 



(1) 



20 24 




20 24 



FIG. 3. The 2D electronic concentrations in occupied sub- 
bands. Solid lines represent theoretical results, marked points 
denote data derived from magnetoresistance measurements. 
Due to strongly field-dependent cyclotron effective masses, 
measurements at different groups of angles were appropriate 
for different subbands and different in-plane field regions. Na 
and Ni, denote experimental results obtained from the Fourier 
transformation of Shubnikov-de Hass oscillations measured at 
low perpendicular magnetic fields {B < 1). 



where m* is the effective mass in the GaAs conduction 
band and the vector potential for applied along the 
y-direction takes a form A = {B\^z, 0,0). The confin- 
ing potential Vconf{z) is constructed using the nominal 
growth parameters and includes Hartree and exchange- 
correlation potentials generated by the free carriers in 
the quantum wells. The Hartree potential is derived 
from the 2;-dependent density of electrons, by nu- 

merical solution of the Poisson equation. The exchange- 
correlation terrtL is calculated within the local-density 
approximation.Ej In each loop of the self-consistent pro- 
cedure we solve numerically the Schrodinger equation 
with the Hamiltonian (|^) to get p{z). Then, a new 
Vconfiz) is constructed, which enters the next loop of the 
procedure until the self-consistency condition has been 
achieved. The resulting energy spectra have the form 



(2) 
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where Ei^x{kx) depends on the sample geometry and on 
the magnitude of the in-plane magnetic field. An index 
i = a,b distinguishes the antibonding and bonding sub- 
band states. 

2D subband concentrations Ni are proportional to the 
area enclosed by corresponding Fermi contours. The par- 
tial densities of states (DOS) gi and the cyclotron effec- 
tive masses m^i are related to the shape of the Fermi 
contours by the following expression: 



1 

2n^ 



dk 



(3) 



The total density of states g = ga + gt- The distortion 
of Fermi lines can be probed experimentally by adding a 
weak perpendicular component of the magnetic field. For 
B± <^ , the quantization of the in-plane component of 
the electron motion can be described in terms of quasi- 
classical Landau levels with the cyclotron effective mass 
given by Eq. (|). The de generacy of the spin-unresolved 
levels is 2\e\B±/h. 

The evolution of Fermi surfaces in is illustrated in 
Fig. 1^ for sample B. Fig. ^(a) presents theoretical Fermi 
contours calculated for several selected values of , cor- 
responding cyclotron effective masses are shown in Fig. 
^(b). The deviations from the zero-field circles reflect 
the field-induced changes of Ei^^ikx) described for the 
case of the simple tight-binding model in Section |[ The 
antibonding Ea^xikx) is an asymmetric function of kx, 
narrower than the free-electron parabola, with the mini- 
mum (bottom) moving to the Fermi energy which crosses 
at _B|| = -Bc,i. The corresponding Fermi contour acquires 
the shape of a "lens" , i.e. the oval with the longer axis 
oriented in the /cj^-direction. Due to this type of defor- 
mation the antibonding cyclotron effective mass m^a is 
a decreasing function of B\\ . 

The field-induced local maximum developing in the kx- 
dependence of Eb_xikx) causes elongation of the bonding 
Fermi contour in the fcj;-direction. For higher fields, the 
Fermi line evolves into an asymmetric "peanut" shape. 
A neck connecting left and right parts of the "peanut" 
breaks as the saddle point reaches the Fermi energy at 
B\\ = Bc.2- For _B|| < Bc,2 the bonding effective mass 
nic.b grows with and diverges at the saddle point. 
Above i?c,2 7 the Fermi contour splits into two approxi- 
mately elliptic lines. Very similar results are obtained 
for sample B. 

The position of an electron described by the state 
\i,kx,ky) in one or the other well is given by {z)i^x — 
{i,kx,ky\z\i,kx,ky). The coordinate {z)i^x is related to 
the fca;-component of the wave vector by 
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FIG. 4. (a) Theoretical evolution of Fermi contours in the 
in-plane magnetic field, (b) Self-consistently calculated cy- 
clotron efi'ective masses. 



where to — |e|-B||/m* and {v)i^x is the x-component of the 
in-plane group velocity. These relations imply that the 
larger, almost circular Fermi contour describes the elec- 
tron layer in the left rectangular well while the smaller 
elongated oval corresponds to states in the right trian- 
gular well. The cyclotron effective masses in two wells 
are different: above Bc.2 the mass for the left well drops 
quickly close to the zero-field value m*, the mass for 
the triangular well also drops, but to substantially larger 
value. 

Both the partial densities of electrons Paiz), Pb{z) and 
their sum p{z) — Pa{z) + pb{z) are influenced by in-plane 
magnetic fields. In our structure with a hard-wall bar- 
rier the densities have sharp minima in the barrier and 
peaks inside the wells. The position of peaks is almost the 
same for all densities and practically field-independent, 
the magnetic field effect is represented by changes of rel- 
ative magnitudes of the peaks. We will characterize the 
transfer of electrons between wells by Zi{B^^), the field- 
dependent centroids of electron densities. Results calcu- 
lated for sample B are shown in Fig. ^. The center of mass 
of p(z) is in the middle of the structure at -B|| = and 
moves closer to the interface as the field increases. This 
behavior is typical for single-junction structures Below 
i?c,i) the centroid Za of antibonding electrons is inside the 
triangular well for all fields. At zero-field, Za is far from 
the barrier. The growing Bu transfers antibonding elec- 
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trons from the triangular to rectangular well and, conse- 
quently, Za is shifted towards the barrier. The electrons 
from the bonding subband exhibit opposite behavior in 
this range of fields: the magnetic field empties the an- 
tibonding subband and corresponding electrons, mostly 
from the triangular well, become a part of the bonding 
subband. As result, Zb moves from position close to the 
interface to the position near the barrier as increases. 
Above Bc.i, when only the bonding subband is occupied, 
the transfer of electrons from the triangular to rectangu- 
lar well continues. This is the case also above Bc^2 when 
the bonding Fermi contour is splitted into the left and 
right parts, as seen in Figs. |^ and ^. 
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FIG. 5. Theoretical centroids Zi, i = a,b, of densities Pa{z), 
pb{z) and their sum p{z) — pa{z) + Pb{z) as functions of 
in-plane magnetic field. Note the different scales for the pos- 
itive and negative parts of the vertical axes. 

The characteristic of the samples A and B derived from 
the magnetoresitance experiments are collected in Figs. |l| 
and ^, together with results of the self-consistent numer- 
ical simulation. In the case when theory cannot be di- 
rectly compared with experiments only the theoretical 
results for sample B are presented. 

We associate the sharp minima in the experimental 
magnetoresistance traces, presented in Fig. with the 
sudden decrease of the density of states at i?c,i that con- 
tributes to reduction of the scattering rate at the Fermi 
energy. The difference between the field-dependence of 
cyclotron effective masses rric a and TOc,b, calculated for 
sample B and shown in Fig. |j(b), explains why the SdH 
oscillations of electrons in the antibonding subband are 
seen at lower tilt-angles than the oscillations of electrons 
from the bonding subband. 

The increase of theoretical m^b at in-plane fields be- 
low i3c,2 is related to the strong distortion of the Fermi 
contour, as seen in Fig ^. Only the bonding subband is 
occupied for i?c,i < -B|| < -Be, 2 and, therefore, the SdH 
oscillations are periodic in 1 / B± . The decrease of their 
amplitude for _B|| approaching i?c,2 is an experimental 
confirmation of calculated sharp increase of TOc,fc- The 



theoretical value for the critical in-plane field precisely 
matches with i?c,2 obtained from experiment. The peak 
in the measured in-plane field magnetoresistance at Bc.2 
can be understood as a consequence of the divergence of 
the DOS at Fermi energy or, in other words, of a zero 
group velocity of electrons at the saddle point. 

The experimental magnetoresistance oscillation is 
clearly associated with van Hove singularities in the 
dependent DOS. However, apart from these strong fea- 
tures Rxx does not follow g for in-plane fields less than 
Bc,2 (see Fig. Q). The differences between AR^x/Rxx and 
Ag/g curves stem from the in-plane field suppression of 
the coupling between wells, the transfer of electrons from 
the triangular to rectangular well, and the changes in the 
nature of electronic scattering. 

At _B|| > -Be, 2 7 the Fermi contour splits into two ap- 
proximately elliptic lines. The larger, almost circular 
Fermi contour is for the rectangular well while the smaller 
elongated oval corresponds to states in the triangular 
well. Since electrons are localized in one or the other 
quantum well, the scattering rate becomes nearly inde- 
pendent of and ARxx/ Rxx oc Ag/g applies, as seen 
from Fig. |l|. Large cyclotron effective mass and low 
concentration of electrons in the triangular well explain 
while these states are not detectable by the SdH measure- 
ment; at most two weak oscillations can occur in avail- 
able magnetic fields. The calculated density of electrons 
in the rectangular well is again in excellent quantitative 
agreement with SdH data. 

IV. CONCLUSIONS 

We have performed an experimental and theoretical 
study of asymmetric GaAs/Alo.sGao.yAs heterostructure 
consisting of coupled rectangular and triangular quantum 
wells. The positions of a strong minimum and maxi- 
mum in the measured in-plane field dependent resistance 
of the 2D electron system match precisely with calcu- 
lated van Hove singularities in the DOS. Theoretical pre- 
dictions for -dependent occupations of electronic sub- 
bands and quantum wells are in quantitative agreement 
with SdH data recorded at magnetic fields slightly tilted 
from the 2D layer plane. We conclude that the numerical 
self-consistent-field technique provides a realistic descrip- 
tion of orbital effects of the in-plane field on 2D electron 
systems confined in semiconductor heterostructures with 
general growth-direction geometry. 

To achieve similar level of accuracy for valence- 
band states with strong spin-orbit coupling represents 
a formidable challenge for future theoretical work. Par- 
ticularly intriguing is the derivation of local-spin-density 
approximation for many-body states with unequal pop- 
ulation of different spin-subbands. Recent remarkable 
studies of the metal-insulator transition in high mobility 
GaAs hole systems certainly give a strong motivation for 
pursuing this kind of research. 
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